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It's On: Early Interpretations of ATLAS Results in Jets and Missing Energy Searches 
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The first search for supersymmetry from ATLAS with 70 nb -1 of integrated luminosity extends 
the Tevatron' s reach for colored particles that decay into jets plus missing transverse energy. For 
gluinos that decay directly or through a one step cascade into the LSP and two jets, the mass range 
vrtg < 205 GeV is disfavored by the ATLAS searches, regardless of the mass of the LSP. In some cases 
the coverage extends up to wig c± 295 GeV, already surpassing the Tevatron's reach for compressed 
supersymmetry spectra. 



I. INTRODUCTION 

Jets plus missing energy is one of the most promising 
channels in the search for physics beyond the Standard 
Model at hadron colliders. With the preliminary results 
from ATLAS on jets plus missing energy searches for su- 
persymmetry PQ, new limits on color octet particles that 
subsequently decay into a pair of jets and a long lived 
neutral particle can be set. Although only 70 ±8 nb -1 of 
integrated luminosity has been analyzed, that is already 
sufficient to infer new constraints on models with light 
colored particles. 

The most robust bounds on new color octet particles 
come from LEP2's measurements of thrust. These mea- 
surements place a model independent lower bound on a 
new color octet particle of 51 GeV [2 . Above this en- 
ergy, the limits are model dependent and derive from 
searches that trigger on jets and missing energy. In fact, 
most recent searches have focused on models closely re- 
lated to mSUGRA [3 . Within the many assumptions of 
mSUGRA lies the condition of gaugino mass unification. 
This requirement imposes that throughout all of the ac- 
cessible mSUGRA parameter space, the LSP is a bino, 
the fermionic partner of the hypercharge gauge boson, 
and the gluino is approximately six times heavier than 
the bino (see [4] for recent Tevatron searches). There- 
fore, exclusion limits that are stated in the mSUGRA 
parameter space never effectively vary the kinematics of 
the signal. 1 

Even within supersymmetric models with gauge cou- 
pling unification, gaugino mass unification is not an auto- 
matic outcome. For instance, in mirage mediation which 
mixes gravity mediation with anomaly mediation, the 
gauginos can be degenerate at the TeV scale [BJ. In non 
supersymmetric models that also have a gluino-like ob- 
ject in their spectrum, such as Universal Extra Dimension 
models (UEDs) with the lightest KK gluon, gaugino mass 
unification is typically not satisfied [7 J. Therefore, most 
jets plus missing energy studies do not directly apply to 
these theories. 



UA2 provides some of the most extensive discussion on gluino 
limits away from gaugino mass unification 



Moreover, in mSUGRA parameterizations the gluino 
and squarks' masses are related, and analyses based on 
such parameterizations rely on the production of squarks 
as well as gluinos to place indirect bounds on the gluino 
mass. 

Model independent studies reinterpreted the Teva- 
tron's searches for jets and missing energy in a more 
general setting [8]. Models where the gluino decays 
into jets plus the LSP have a universal lower bound 
of m~g > 130 GeV (stronger constraints apply to spectra 
where the gluino and LSP have hierarchical mass split- 
tings). This universal limit is far below the normally 
quoted 400 GeV limit on gluino-like objects. Given such 
modest current bounds - a 51 GeV hard lower limit on 
the gluino mass and an approximate 130 GeV lower limit 
from dedicated searches - the integrated luminosity ac- 
cumulated at the LHC in its early running is already 
sufficient to extend the present reach for light colored 
octets. 

Low luminosity searches allow modest cuts to be placed 
on missing transverse energy, helping the searches' effi- 
ciencies for discovering compressed spectra. ATLAS per- 
formed four separate searches for excesses in jets plus 
missing energy: monojets, dijets plus missing energy, 
three jets plus missing energy and four jets plus miss- 
ing energy. These searches required a leading jet with 
Pt > 70 GeV and $t > 40 GeV. This low $t cut is par- 
ticularly sensitive to theories with somewhat compressed 
spectra and improves over the Tevatron's limits. 



II. MODEL AND SIMULATION 

This article considers a limit of the MSSM where the 
gluino is the only accessible colored particle, and decays 
via an off-shell squark into two jets plus the lightest su- 
persymmetric particle (LSP). More specifically, the the- 
ory consists of a color octet Majorana fermion, g, with no 
SU(2)lxU(1)y charges, and a lighter Majorana fermion 
with no Standard Model gauge quantum numbers, x° • 
X° is stable on detector timescales and may be the dark 
matter of the Universe, g can decay into x° plus two 
jets or to x° an d a gluon. For simplicity, it will be as- 
sumed that the branching ratio into light quarks is unity, 
Br(g — » x°QQ) = 1- The decay is mediated by a dimen- 
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C decay = ^( X °qa)T Aa b (Q b ~g A ), 



(1) 



which can be generated by integrating out a color triplet 
scalar. The lifetime of g is approximately 



{m~ g -m x of 
4ttA 4 



(2) 



which leads to a prompt decay so long as A < 10 TeV 



for 



m x o + 10 GeV. There is no a priori relation 



between the masses of x° or 9- X° mav be very light 
without any constraints arising from LEP, and the only 
model independent constraint on g is that it should be 
heavier than 51 GeV. 

Models with approximate g~x° mass degeneracy are 
particularly challenging for jet plus $t searches. In this 
case, the decay products of g are soft and buried in back- 
ground. In the degenerate limit, the most efficient way 
to detect g production is by looking for radiation of addi- 
tional jets. At the Tevatron, pair produced g's plus radia- 
tion gives rise to events with low multiplicity jets plus $t- 
In particular, monojet searches can be effective at dis- 
covering these topologies [8]. However, monojet searches 
are typically exclusive and place poor bounds away from 
the degenerate limit. For instance, CDF places a sec- 
ond jet veto of Et j 2 < 60 GeV and a third jet veto of 
Et j 3 < 20 GeV [21]. As the mass difference between g 
and x° increases, the efficiency of such cuts diminishes. 
In the non-degenerate limit, the most suitable searches 
have higher jet multiplicity. However, the cuts applied 
on the monojet and multijet searches performed so far 
were sufficiently strong that they left a gap in the cov- 
erage of the intermediate mass-splitting region [8 j. The 
present bound on nig only extends above 130 GeV for 
m x o < 100 GeV. The LHC cross section for gluinos just 
above this limit is of the order of a few nanobarns. There- 
fore, limits can be improved by the LHC with remarkably 
low luminosity and early discovery is potentially achiev- 
able. Unfortunately, no excesses were observed in pQ, so 
only new limits could be inferred. 

In this work, the efficiencies of the cuts applied by 
ATLAS' recent search are extracted through a Monte 
Carlo study. These efficiencies depend on nig and m % o 
and are necessary to calculate limits. The signal is 
calculated using MadGraph 4.4.32 [9 , matching parton 
shower (PS) to additional radiation generated through 
matrix elements (ME) using the shower k± matching 
prescription from [10J, which closely mimics the MLM 
matching scheme [11 . The shower k± scheme allows for 
the matching scale, Q C ut? to be set equal to the matrix 
element cutoff scale, Q^f , and it more directly samples 
the Sudakov form factor used in the parton shower. In 
the region where g and x° are nearly-degenerate, the ad- 
ditional radiation is crucial in determining the shape of 
the $t distribution and hence how efficiently the signal 
is found. A matching scale of Qcut = 100 GeV is adopted 
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FIG. 1: pr spectrum for the 1st, 2nd and 3rd hardest jets 
in the events (blue, red and green, respectively) when ISR 
matching is included (solid lines) vs when all radiation is 
accounted for by the parton shower (dashed lines). An in- 
tegrated luminosity of C — 70 nb _1 is assumed, as well as 



masses ms — 100 GeV and ra v 



90 GeV. 



for the signal and the matrix elements for the following 
subpro cesses are generated: 2g + 0j, 2g + lj and 2g + 2 + j. 
When performing MLM matching all higher multiplicity 
jet events are generated through parton showering. The 
choice of a s (/i) in MLM matching is not varied. Fu- 
ture studies should document the effect of scale choice 
in matching and incorporate this theoretical uncertainty 
into limits. 

The importance of matching in describing hard radi- 
ation is particularly pronounced for models with nearly 
degenerate spectra. That is illustrated in Figure [I] for 
a benchmark mass point vrtg = 100 GeV and m % o 
90 GeV, which contrasts the pr spectrum of the 3 hard- 
est jets in the events when matching is included vs when 
all radiation is due to the parton shower only. 

The parton showering is performed in PYTHIA 6 . 4 [12] . 
PYTHIA also decays g —> qqx° , hadronizes the events and 
produces the final exclusive events. These events are then 
clustered using a cone-jet algorithm with R = 0.7 with 
PGS4 [13 which also performs elementary fiducial vol- 
ume cuts and modestly smears the jet energy using the 
ATLAS-detector card. 

ATLAS' search does not use proper $t in their analy- 
sis, instead it uses missing energy at the electromagnetic 
scale, ^tem- The relation between $t and ^tem is 
shown in Fig. 7 of [14]. In this paper, it is approximated 
as 



$T EM 



1.5-# T /2100 GeV 



(3) 



where Ht is the sum of the transverse energies of the 
jets in the event. This effectively raises the $t cut to 
approximately 50 GeV. 

In order to validate this modeling of ^tem, the SU4 
mSUGRA model shown in [1 is reproduced. The SUSY 
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Cut 


Topology 


1 j + ^ T 


2+ j + $ T 


3 + j + $t 


4+ j + $ T 


1 


PTl 


> 70 GeV 


> 70 GeV 


> 70 GeV 


> 70 GeV 


2 


PTn 


< 30 GeV 


> 30GeV(n = 2) 


> 30GeV(n = 2,3) 


>30 GeV(n = 2-4) 


3 


-^TEM 


> 40 GeV 


> 40 GeV 


> 40 GeV 


> 40 GeV 


4 


PT£ 


< 10 GeV 


< 10 GeV 


< 10 GeV 


< 10 GeV 


5 


Ac/)(j ni ^tem) 


none 


[> 0.2, > 0.2] 


[> 0.2, > 0.2, > 0.2] 


[> 0.2, > 0.2, > 0.2, none] 


6 


^TEM/M e ff 


none 


> 0.3 


> 0.25 


> 0.2 




^Vp r ed 


46t£ 


6.6 ±3.0 


1.9 ±0.9 


1.0 ±0.6 




N 0hs 


73 


4 





1 




e x a(pp -> ggX)\ 95% c .l. 


663 pb 


46.4 pb 


20.0 pb 


56.9 pb 



TABLE I: Searches in [T used to set limits in this article. The 95% C.L. on the production cross section times efficiency of the 
cuts, e x a(pp — »• ggX), follow from folding in the uncertainties in the luminosity and background. 



Les Houches Accord parameter card [15] for SU4 is cre- 
ated with a spectrum calculated with SuSpect 2 .41 [16] 
which matches the spectrum reported in [17] to 5% accu- 
racy. The decay card for SU4 is calculated with SDECAY 
[18], interfaced with SUSY-HIT [19], and finally the cross 
sections are generated with MadGraph and decayed, show- 
ered and hadronized in PYTHIA. The total SUSY produc- 
tion cross section is normalized to the NLO value used 
in pQ in order to compare efficiencies and shapes of dis- 
tributions. 

The $t em distribution for the SU4 benchmark model 
calculated with PGS4 is compared against the $t em dis- 
tribution in pQ . Using the $t em parametrization in Eq. 
[3j the normalization of the PGS4 distribution is lower by 
7% than the corresponding values in Fig. 10a of pQ. The 
discrepancy in the individual bins of the $t em distribu- 
tions can be quantified by the RMS variance, 
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0.1, (4) 



where iVbms is the number of bins, and nf^ TLAS (nf GS4 ) 
are the number of events in the i th bin of the ATLAS 
(PGS4) distribution. The absolute efficiency of the fixed 
$t em cut (Cut 3 of Table [l| in PGS4 differs from the 
ATLAS' absolute efficiency by ±12%. 

Fig. 8 of [1 shows the fractional $t distribution for 
SU4 (used in Cut 6 Table [j]). The fractional $ T is defined 
as $T/M e ff, where M e ^ is given by 



(5) 



Here n is the multiplicity of the channel, e.g. n = 3 for 
the 3 + j+^t channel. The difference in the absolute nor- 
malization of these distributions is 7% between ATLAS 
and PGS4. The $t em parameterization in Eq. [3] tends 
to underestimate the SU4 fractional $t distribution, but 
reproduces the efficiency of the fractional $t cut to 15% 
accuracy. The fractional $t cut is more sensitive to the 



modeling of ]£t em and depends upon the exact kinemat- 
ics of the signal. Therefore, ultimately the results of this 
paper should be validated by the ATLAS collaboration. 

The same procedure used for the calibration process 
is applied to the study of the models with new light col- 
ored particles considered here. After the simulated signal 
events are processed with PGS4, the cut efficiencies for the 
PP ~ > 9§X processes are inferred. When combined with 
the results from the ATLAS searches, a limit can be set 
on the maximum production cross section for a given vrtg 
and ra£. 

The cuts applied by ATLAS' early supersymmetry 
searches pQ are listed in Table [l| There are several uncer- 
tainties folded into the limits. The first is the measured 
luminosity, C = 70 ± 8 nb _1 . For the analysis in this 
paper it is assumed that the luminosity is described by 
a Gaussian distribution. The second uncertainty is the 
background expectation. Ref. [1 used background esti- 
mates normalized at low $t to calculate the QCD con- 
tribution to the high $t signal regions. Due to the large 
estimated errors, the number of expected background 
events is assumed to follow a log- normal distribution. 
There is a weak correlation between the background and 
luminosity errors but for simplicity these uncertainties 
are taken to be independent. Both of these uncertainties 
are convolved with a Poisson distribution for signal plus 
background to derive the 95% confidence limits (C.L.) 
on the cross section times efficiency and are also listed in 
Table H 

Since the signal events typically contain four jets, 
higher multiplicity searches are expected to provide the 
best limits. But in practice, the 3 + j + ~$t had the best 
exclusion power, and the 4 + j + $t channel was the least 
sensitive among the multijet searches by factors of several 
because the Prj 4 > 30 GeV cut is not efficient. For in- 
stance, for a spectrum with rrig ~ 300 GeV and m x o ~ 0, 
the typical energy of each of the four jets is ~100 GeV. 
However, since the jets are approximately isotropically 
distributed, the chance of getting all four jets with p T 
above 30% of their momentum is only half as likely as 
having only three out of the four jets passing this cut. 
Moreover, no events were observed in the 3 + j+^t chan- 
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nel versus one event in the 4 + j + jS? channel. 

Near the degenerate limit, the 2 + j + JE? search was 
the most effective and set the tightest limit. The lj + ]£t 
search did not provide competitive limits because of the 
large number of observed events. 

The limits placed on the production cross section al- 
low one to infer the excluded parameter space of specific 
models. The most straight-forward of those is the MSSM 
limit where the gluino and the LSP are the only currently 
kinematically accessible particles. Fig. [2] shows the next- 
to-leading order (NLO) cross section for pp — > ggX calcu- 
lated with Prospino 2 . [20 as a function of the gluino 
mass. 
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FIG. 2: The NLO QCD cross section for pp — > ggX as a func- 
tion of the gluino mass. A cross section of a — 1/(70 nb -1 ) 
corresponds to rrig — 395 GeV and a — 10/(70 nb -1 ) corre- 
sponds to wig = 265 GeV. 



The combined results from all the channels are shown 
in Fig. [3] where only the most restrictive limit on the 
the production cross section is used. Fig. [3] displays the 
contours of the maximum allowed production cross sec- 
tion in the nrig-m x o parameter space. There is a slight 
loss of sensitivity as m x o — » while holding rrig fixed be- 
cause the cut on the fractional missing energy has lower 
efficiencies for larger mass splittings. This is because as 
X° becomes lighter, JE? grows, but the sum of the jet 
energies grows faster, thereby decreasing the fractional 
missing energy. 

The searches performed by ATLAS also allow new 
bounds to be placed on pair produced color octets that 
decay through a one step cascade into jets plus JSt- How- 
ever, the limits have reduced sensitivity because cascade 
decays reduce the ]£t in the events [22]. Consider the 
following simplified model as a way of studying the ef- 
fects of cascade decays on the search for color octets. As 
before, there is a color octet Majorana fermion, g. In 
addition to a neutral, Majorana fermion x°, there is an 
intermediate charged fermion, x ± - The x ± decays to 
X°W ± (*\ where the W ± is either on-shell or off-shell de- 
pending on whether the mass splitting between x ± an d 
X° is greater than or less than m w ± . g decays with 100% 
branching ratio to q / qx ± , with equal probabilities for x + 
and x~ - The mass of x ± is taken to be the arithmetic 
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FIG. 3: 95% C.L contours of the maximum allowed produc- 
tion cross section a(pp — ► ggX) in the rrig — m x o mass plane, 
for g directly decaying to x°33- The contour values are speci- 
fied in the right color scale. The dark line corresponds to the 
exclusion boundary for models where the gluino is produced 
through QCD alone with an NLO cross section (i.e., with 
all squarks decoupled so that there are no t-channel squark 
exchange diagrams). The dashed-lines delimit the excluded 
parameter space of different models where a(pp — ► ggX) is 
given by a simple rescaling of the NLO-QCD cross section. 
The red line is the current estimate of Tevatron limits taken 
from 8 . The blue line denotes a sample mSUGRA spectrum 
where g is the gluino and x° is the bino. 



average of rrig and m % o, m x ± = (rrig + m % o)/2. This 
choice of spectrum is conservative because it reduces the 
$t nearly maximally. Fig. [4] shows the limits on this sim- 
plified model. When the W ± goes on-shell, the energy 
put into invisible states decreases and there is a precipi- 
tous loss of sensitivity. The ATLAS search improves the 
limits in this case as well, placing a universal bound of 
rrig ~ 205 GeV on the gluino mass. 



III. DISCUSSION 

Based on ATLAS' recent supersymmetry search, new 
limits for a color octet particle that decays into jets and 
a missing neutral particle were calculated. Surprisingly, 
with remarkably low luminosity these searches surpass 
the previous limits from Tevatron's analyses. This is 
possible due to a combination of the low cuts on missing 
energy, the quality of the detectors and the level of un- 
derstanding that the experiments have already gained on 
jet energy calibration. In the future, 6-jets plus $t will 
be another important subclass of topologies to search for 
and may quickly improve limits over the Tevatron. 

This study of the 70 nb -1 ATLAS' search applies to 
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m- g (GeV) 

FIG. 4: 95% C.L. contours of the maximum allowed pro- 
duction cross section a(pp — ► ggX) in the vrtg — m x o mass 
plane, for a one-step cascade decaying gluino, g — > q'qx^ ~~ 
q / qX°W ±< "*\ The contour values are specified in the right 
color scale. The purple line corresponds to m x ± = m x o + 
m w ± ; below it the W ± goes on-shell. The dark line cor- 
responds to the exclusion boundary for models where the 
gluino is produced through QCD alone with an NLO cross 
section (i.e., with all squarks decoupled so that there are no 
t-channel squark exchange diagrams). The dashed- lines de- 
limit the excluded parameter space of different models where 
a(pp — > ggX) is given by a simple rescaling of the NLO-QCD 
cross section. The red line is the current estimate of Teva- 
tron limits taken from [§]. The blue line denotes a sample 
mSUGRA spectrum where g is the gluino and x° is the bino. 



many supersymmetric models with a light gluino in the 
spectrum (see [23 for a more general discussion of limits 
on non-mSUGRA-like supersymmetric theories). It also 
applies to many non-supersymmetric models with a dark 
matter particle such as theories of five dimensional UEDs, 
six-dimensional UEDs with scalar adjoints [24 , and little 



Higgs models with additional colored composites [25] . 

While color octets that decay into jets plus missing 
energy are a common feature in models of new physics, 
there are numerous nearby models that also merit study. 
These include color triplets that decay into a single jet 
and missing energy, or new colored particles that cas- 
cade decay down to the dark matter. The results pre- 
sented here are applicable to these alternative scenarios 
as well; in fact, they should apply to any theory where 
pair-produced new particles each decay to two jets plus a 
long lived neutral particle, so long as the spectrum is not 
approximately degenerate. In the degenerate limit, ad- 
ditional QCD radiation in the events becomes important 
and therefore the color of the new particle is relevant. 

The LHC is a huge leap forward and even early anal- 
yses can provide coverage of regions that have not been 
covered by the Tevatron or earlier experiments. Ulti- 
mately the exact results in this article need to be re- 
produced and refined by ATLAS to account for detector 
effects. CMS also has great promise in this channel and 
their particle flow algorithms give spectacular measure- 
ments of jet energy. For instance, even CMS' calibration 
data [26] is able to slightly extend the Tevatron's lim- 
its. Going forward, each experiment has accumulated 
roughly 1 pb _1 of integrated luminosity and therefore 
there are great possibilities for discoveries in 2010 as the 
LHC's reach rapidly expands into new kinematic regimes. 
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